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We propose a novel symmetrization procedure to beat decoherence for oscillator-assisted quan- 
tum gate operations. The enacted symmetry is related to the global geometric features of qubits 
transformation based on ancillary oscillator modes, e.g. phonons in an ion-trap system. It is shown 
that the devised multi-circuit symmetrized evolution endows the system with a two-fold resilience 
against decoherence: insensitivity to thermal fluctuations and quantum dissipation. 
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PACS numbers: 03.67.Lx, 03.65.Yz, 03.65.Fd 

A key requirement to implement quantum computa- 
tion consists in the ability to perform, with high accu- 
racy, some set of universal gate operations for quantum 
bits (qubits) Q. In the real world, a quantum system 
is never isolated, and the decoherence due to the sur- 
rounding environment is regarded as the main obstacle 
to practical realization of quantum information process- 
ing. A variety of techniques, based on the underlying 
mechanism of symmetry, have been proposed to protect 
quantum information, including quantum error correct- 
ing codes 0, decoherence-free subspaces bang-bang 
dynamical decoupling 0, and noiseless subsystems [EE!- 

Most of these strategies rely, in a way or the other, 
on symmetry properties of the dynamics of the plain in- 
teractions of quantum information manipulation [y] . For 
instance, in the proposed schemes for dynamical decou- 
pling Q , the qubit-qubit coupling is often invariant with 
respect to the symmetrization procedure exploited for 
the elimination of the qubit-environment couplings. On 
the other hand, it has been recently argued that robust 
quantum computing might be also achieved by a class 
of geometric transformations 7]. In particular, geomet- 
ric quantum manipulation seems to possess a certain in- 
triguing resilience to classical noises, i.e. stochastic fluc- 
tuations of the control parameters, due to its inherent 
global feature of the operation Q . 

In this paper, we shall contrive a novel symmetriza- 
tion procedure to beat decoherence for a sort of geomet- 
ric quantum gates based on the use of ancillary oscillator 
modes. For the typical ion-trap system, the global gate 
operation assisted by a phonon mode and the related 
time-reversal symmetry will be exploited. At variance 
with the previous schemes, the symmetry involved in our 
proposal will not directly show up at the level of dynami- 
cal interactions, i.e. Hamiltonian operator, but will relate 
to the geometric nature of the overall, finite-time evolu- 
tion, i.e. unitary operator. The use of this sort of global 
transformations to implement quantum gates, and their 
extra advantages, have been already discussed by several 
authors [Io| . However, the symmetry associated with 
these geometric operators and its potential relevance to 
the struggle against decoherence was never revealed. 

Let us focus, for the sake of concreteness, on quantum 



information processing based on ion-trap systems. Quan- 
tum manipulation on this kind systems constitutes one of 
the most promising approaches to quantum computation 
and it has been intensively studied during the past years 
[ToL ITU IT^ . In detail, qubits in the system are stored 
in internal atomic states controlled by laser beams and a 
controlled interaction between them is obtained via the 
coupling with ion vibrational states. Although the cou- 
pling of the qubits with the phononic degrees of freedom 
will generally lead to dependence of the gate operation 
on the motional states, recent studies show that such a 
drawback can be avoided by resorting to a control pro- 
cess enacting a gate with global features . The main 
idea of the scheme can be illustrated as below. Suppose 
that, by appropriately choosing the laser intensities and 
detunings, the following coupling of two ion qubits with 
the oscillator degree of freedom can be obtained 



H(t) = [a(t)a + a* {t)oft]J y , (1) 

where the collective Pauli operator J y = ay + <jy ac- 
counts for the internal states of the ion qubits and a(t) 
characterizes the interaction of the qubits with the vibra- 
tional state. When a(t) has the form of certain periodic 
function such that in the time interval J Q a(t)dt = 0, 
the corresponding evolution generates a global transfor- 
mation, e ~ lA ^ J y with A(T) = ifa a(t)f*a*(t')dt'dt, 
which acts only on the internal states of the ions. The 
simplest case is given by the periodic step pulses, which 
leads to the following gate transformation 



U C (T) = PI{t)PI(t)P1(t)P p _{t) = e 



-i2af,T z Jf, 



(2) 



where ao and t = jT stand for the pulse amplitude and 
length respectively; P±(t) — e ±l ^ a ° TxJ y and P±(t) — 
e ±iV2a Tpj y (ies^-ipg ^e unitary transformations for cor- 
responding pulse interactions. [We have used the nota- 
tion x = -^=(a + eft) and p = -4= {aft — a)]. Note that 

the transformation (J2J contains no operators with a non- 
trivial action over the vibrational degrees of freedom. 
The advantage of this gate operation is evident: it is in- 
dependent of the vibrational states, hence is insensitive 
to the ion temperature. 
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FIG. 1: The path traversed in phase space of the circuit evo- 
lution: (a) the case of step pulse interaction; (b) the circular 
case with a(t) = ao cos |^ + iao sin 7p. 



One more feature of the above described operation, 
which will be at the very basis of the scheme to beat de- 
coherence proposed in this paper, is the invariance of the 
global action under the time reversal of the interaction: 
Hc{t) —> H(j{t) = —Hc(t). This property can be easily 
checked in the case of the step pulse interaction. Indeed 
one has: U (T) = P*(t)PP(t)P*{t)PI(t) = U C {T). 
Represented in Fig. 1 are the corresponding circuits 
in the oscillator phase space, including the above step 
pulse interaction and also the circular case with a(t) = 
ao cos ^ + iao sin . The invariance of the global action 
under time reversal has, of course, a simple geometric in- 
terpretation: the corresponding paths traversed in the 
phase space, being related by the map a — > —a, indicate 
an anti-symmetric track geometry. 

In real-life systems, the qubits and the oscillator de- 
grees of freedom cannot be perfectly isolated. Here we 
shall focus on the case in which decoherence effects will 
arise by a thermal dissipation of the oscillator mode only. 
This latter will spoil the global action of and induce 
qubit-oscillator entanglement which amounts to decoher- 
ence. The relaxation due to the coupling between the 
oscillator and the environment will be described by the 
following master equation 



d t p = -i[H(t),p] + ^n(2apa) - {a ! a,p}) 
:= (C + C t h)p, 



(3) 



where K denotes the relaxation rate. The Liouville op- 
erator contains two terms C — Co + Cth with Cop — 
—i[H(t), p\. Let us first examine, for transparency, the in- 
fluence of dissipation on the no-jump trajectory, namely, 
the evolution generated by the following non-Hermitian 
Hamiltonian 



H(t) = H(t) 



(4) 



Let us consider the circuit of step pulses applied to the 
system. If the relaxation rate k is small, i.e. k <C 
a Q ~ 1/t, then the non- unitary transformation gen- 
erated by every pulse, P±' p (t), can be approximately 
expanded by resorting to the Baker-Hausdorf relation 



way it is not difficult to realize that, to the first order in 
rt/aoj the sequence of four step pulses gives rise to the 
following transformation 



U C (T) = Pf(r)P^(r)Pf(r)P^(r) 



(5) 



g — i2a 2 T 2 Jy e ~ V2ku t 2 xJ v gT/2K,a T 2 pJ y e ~ 2kto) a 

Clearly, the appearance of the middle two factors in the 
last term of Eq. (JjJJ leads to the entanglement of qubits 
with the oscillator degrees of freedom, hence results in de- 
pendence of the transformation on the vibrational state. 

The key observation is now to notice that the above 
fault can be effectively suppressed by performing a fur- 
ther time-reversed pulse sequence. One can readily ob- 
tain that the no-jump trajectory of the time-reversed 
pulse sequence will generate the transformation 



U e (T) = PI{t)PI{t)P*{t)P*{t) 

^ — i2a 2 T 2 J 2 </2~K,a T 2 x,Jy — \Z2Ka T 2 pJ v 



(6) 



From Eqs. (J5J and ©, it can be seen that the influence 
of dissipation in the two circuits counteract each other. 
Therefore, a gate operation resistant against dissipation 
can be enacted by combining these two opposite pulse 
circuits 



U(T 2 ) = U c (T 2 /2)U c (T 2 /2) 



e -t2( ao rYj 2 e 



(7) 



where T 2 = \[2T is the whole time period of the double- 
loop process. Since the two factors of the last term of 
Eq. (7J act individually on the qubits and the oscillator 
degree of freedom, the dissipation-induced dependence 
of the gate operation on the oscillator state has been 
removed. 

It is important to stress that the cancellation of the 
unwanted influence of dissipation achieved by the above 
procedure is essentially a sort of group-symmetrization 
decoupling process Q. Notably, in the present system, 
the enacted Z2-invariance of the dynamics under the map 
J y — > — J y resides in the geometric feature of the global 
transformation J3J; therefore the devised symmetrized 
(anti-symmetrized) evolution (JJJ endows the scheme with 
a two-fold resilience against the noise sources due to the 
coupling with the oscillator: thermal fluctuations and 
dissipation. 

The analysis carried over up to now was limited to 
the evolution in the no-jump trajectory. To get a full 
picture of the dissipative dynamics, i.e. the resolution 
of the master equation, generally requires to resort to a 
numerical approach. Remarkably, as we shall show in 
the following, the master equation of the present system 
© can be solved analytically by using group-theoretic 
techniques at a super-operatorial level. 

Before presenting the solution, it is convenient to in- 
troduce the following notation ■ Let us define the set 
of left operators 0\ and right operators % R : the left 
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operator action on the density matrix is given, as usual, 
by multiplication on left whereas the action of the right 
operators is defined as 

Or x 0% x • • ■ x O n R l p = P O m x ■ ■ • x O 2 x O l . (8) 

By means of these notations, the operators contained in 
the master equation © can be conveniently described. 
For example, the operator Cth, which describes the pure 
relaxation process, can be denoted by Cth = ^k(2K- — 
Kq + 1), where K = a\a L + a R a R , K- = a L a R 
and K+ = ar L a,R. Moreover, in view of the fact that 
these (super)operators generate a simple su(l,l) alge- 
bra: [K ,K ± ] = ±2K±, [K-,K+] = K , the process of 
pure relaxation is given by e Ctht = e~ K - e^~( 1 ~ K °' > e K - . 
The established formalism lead in a rather straightfor- 
ward fashion to the explicit description of the relaxation 
dynamics for a general oscillator coherent state \(3) |14j . 

e Ctht \{3){[3\ = \(3e~i Kt )(Pe-^ Kt \. (9) 

Suppose now that the finite-time evolution governed 
by the master equation @ is given by a map A, i.e., 

pit) = A(i) P (0). (10) 

Note that a natural picture for the system is provided 
by basis of Jj,-eigenstates \l). The time evolution for 
every element of the operator A(t), defined as A" (t) = 
(l\A(t)\l'), can be obtained: 

d t A ll '(t) = (4 1 ' + C th )A lV (t), (11) 

where 

= -il{aa L + a*a\) + il'{aa R + a*a R ). (12) 

It should be clear that Eq. (fTT]) suggests a linear al- 
gebraic structure, including, besides the aforementioned 
su(l,l) algebra, the both the left and right Heisenberg 
algebras {<2£,a^,l} and {a R , a R ,l}. According to the 
group-theoretic approach developed in Ref. 0, such a 
linear algebraic system is exactly solvable. We summa- 
rize below the solution of Eq. ]1H , details of the resolu- 
tion will be reported elsewhere (14| . 

A ll '(t) = e -^"'(*) e -««;(t)4 e -«^e+(t)+ A n'e-(t)i°i. 

where Aw =1 — 1' and the coefficients are given by 

6fc(t) = e**"* f a(t')e ± ^ Kt 'dt', 
Jo 

A u '(t) = iKll' f \Z+{t')\ 2 dt' 
Jo 

Jo 



The above analytical solution fully characterizes the de- 
coherence effects induced by the dissipation of the oscil- 
lator mode, e.g. vibrational relaxation in an ion trap. 
The integrals in the functions £±(i) might clearly be non 
zero; therefore, the transformation described by equa- 
tion 1(13(1 will be inevitably depend on the oscillator 
state. More specifically, let us suppose that the oscil- 
lator state is prepared in a general coherent state |/3) 
and the whole system is initially prepared in the product 
state p{0) = p m (0) <g> \0)<j3\. By substituting A 11 ' (t) of 
1(13(1 into equation 110(1 and tracing over the oscillator de- 
gree of freedom, one obtains the evolution of the reduced 
density matrix of the qubits 

fy) = e -^\t) e lV\ M t)\- 

xe - ! A 1! ,( W _( t )+rci(t))e-i-^ (0); (15) 

where we have taken the element pf n (t) in the J y repre- 
sentation. 

It should now become evident that the dissipation- 
induced effects can be suppressed by the symmetrization 
procedure. The application of the time-reversed pulse 
circuit a{t) pulse will quench the integrations contained 
in £±(i). Specifically, note that the polynomial expres- 

sionsof£±(i) is given by = e*i Kt ^I {k) (t), 
where jW(t) = J* a(t'){^t') k dt' . For the step pulse se- 
quence indicated by Eq. @, one has 7(°)(T) = and 
nonvanishing terms I^ k \T) oc (niT) k of k > 0; but for 
the time-reversal symmetrized pulse sequence indicated 
by Eq. 0, there are /(°)(T 2 ) = I (1) (T 2 ) = 0, leaving 
only higher order terms of k > 1. In fact, such a sym- 
metrization procedure amounts to a decoupling, in which 
the dissipation effects are removed to the lowest-order in 
the scale n/a^. In principle, it is possible to cancel the 
decoherence effects to an arbitrary high order by an iter- 
ative application of the above symmetrization procedure. 
The /c-order decoupling pulse sequence to implement the 
gate operation J2J) is illustrated below 

U {k) {T n ) = U c UcUcU c UcU c U c Uc, (16) 

(i) 

V v ' 

(2) 

(fc) 

where n — 2 k is the number of the pulse circuits and 
T n = \/nT is the total time length of the above pulse 
sequence. 

It is worthwhile to point out that besides inducing the 
mentioned dependence of the qubit evolution on the oscil- 
lator state, dissipation will also result in a deviation from 
the desired gate operation. This latter effect is indicated 
in the coefficient A 11 (t). Note that this undesired influ- 
ence of dissipation, being embedded within the geometric 
action, cannot be removed by the above symmetrization 
procedure. In Fig. 2 we report in detail the fidelities of 
the gate operation for different applications of the pulse 
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FIG. 2: The fidelity of the gate operation [Eq. (T7J] for var- 
ious applications of step pulse sequences. Curves I, II, and 
III correspond respectively to the symmetrization procedure 
with decoupling order k = 0, 1 and 2, in which the number of 
the pulse circuit is taken as n — 1,2 and 4 accordingly. The 
initial state of the oscillator degree of freedom assumes (a) 
P = 2 and (b) = 5. 

sequences as a function of the ratio k/cvo- Fidelity is 
denned as 

F = (V max \p in {T Tl )\* inax ), (17) 

in which T n = i/mr/ao with n the number of the circuits 

and l^max) = e~ l "s j h|00) is the fully entangled state at- 
tained by the ideal gate operation acting on the qubit 
state 1 00). 



In summary, we have contrived a novel symmetrization 
procedure for oscillator-assisted geometric quantum gates 
to suppress the detrimental effects induced by dissipa- 
tion of the oscillator mode. The existence of the symme- 
try is related to the invariance of this class of geometric 
operations under the finite-time reversal transformation. 
We have devised a time-reversal symmetrization proce- 
dure which purports to eliminate quantum decoherence 
due to the oscillator dissipative relaxation to an arbitrary 
high order. The resulting quantum computation scheme 
possesses a two-fold resilient character against decoher- 
ence: insensitive to the oscillator thermal fluctuations 
and immune to quantum dissipation. Note that a similar 
symmetrization procedure was involved in a former refo- 
cusing scheme for holonomic quantum computation [lfij . 
In comparison, the present work suggests a nonadiabatic 
manner to implement a geometrical gate operation ro- 
bust against both classic and quantum noises. Moreover, 
the resistance against dissipation in the present scheme 
has been demonstrated by analytically solving the mas- 
ter equation by means of a powerful group-theoretical 
method at the superoperatorial level. On the practical 
side, we expect our results to be especially relevant for 
ion-trap quantum computing by geometrical gates. 
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